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Using an exact diagonalization approach we show that one- and two-electron InAs quantum dots 
exhibit avoided crossing in the energy spectra that are induced by the spin-orbit couphng in the 
presence of an in-plane external magnetic field. The width of the avoided crossings depends strongly 
on the orientation of the magnetic field which reveals the intrinsic anisotropy of the spin-orbit 
coupling interactions. We find that for specific orientations of the magnetic field avoided crossings 
vanish. Value of this orientation can be used to extract the ratio of the strength of Rashba and 
Dresselhaus interactions. The spin-orbit anisotropy effects for various geometries and orientations 
of the confinement potential are discussed. Our analysis explains the physics behind the recent 
measurements performed on a gated self-assembled quantum dot [S. Takahashi et al. Phys. Rev. 
Lett. 104, 246801 (2010)]. 
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I. INTRODUCTION 

Over the past decade there has been a growing interest 
in the study of the spin-orbit (SO) interaction in semicon- 
ductor low-dimensional systems motivated by the possi- 
bility of coherent spin manipulationii"— The Hamilto- 
nians describing the SO couphng resulting from the in- 
version asymmetry of the material (Dresselhaus^^ cou- 
pling) or the specific structure of the device (Rashbapi^ 
interaction) are not invariant with respect to the rota- 
tion of the spin or the momentum operators separately, 
and consequently spin-orbit-coupled systems posses in- 
trinsic anisotropic properties. This anisotropy has been 
thoroughly studied for delocalizcd systemsi^i^ In partic- 
ular in transport experiments the dependence of the con- 
ductance of a narrow quantum wire on the direction of 
the external magnetic field can be used to determine 
the reciprocal strengths of the Rashba and Dresselhaus 
couplingsi^ The anisotropy of the spin-orbit interaction 
is translated into anisotropic effective magnetic field^ for 
a moving electron modifying the electron spin state. This 
effective magnetic field can be used to perform rotations 
of spin and thus to construct quantum gates^ or a spin- 
field effect transistori^ii^ Moreover, the spin-orbit cou- 
pling is responsible for anisotropic correctionsii to the 
spin swap in a two-qubit quantum gate^^ because it 
results into the precession of spin-packets tunneling be- 
tween the two quantum dotsJ^ 

For electrons localized in a quantum dot the SO cou- 
pling results in avoided crossings (AC) in the energy 
spectra— and spin relaxationii mediated by phonons 
with a relaxation rate dependent on the orientation of the 
external magnetic field.— The energetic effects of the SO 
interaction are usually weak. Only recently SO-induccd 
AC were experimentally measured on quantum-dots that 
were situated in gated nano-wiresi^^ and in gated self- 
assembled quantum dots— The latter experiment stud- 



ied changes of the width of AC for different orientations 
of the magnetic field which extended the previous studies 
that were focused on a comparison of the spin-splittings 
for vertical and in-plane alignment of magnetic ficl d^^i^^ 
in circularly symmetric confinement potentials. 

In the present work we explain the physics underlying 
the observations of Ref. [2ll. To the best of our knowl- 
edge the present paper explains for the first time the 
oscillatory dependence of the width of AC on the direc- 
tion of the in-plane magnetic field. The latter turns out 
to be the consequence of the influence of the individual 
SO couplings and the anisotropy of the confinement po- 
tential. This conclusion is supported by an exact three- 
dimensional calculation of the energy spectra of one- and 
two-electron spin-orbit-coupled quantum dots. 

We show that for quantum dots with confinement po- 
tential elongated in [100] direction for pure Rashba (pure 
Dresselhaus) coupling the AC disappears when the mag- 
netic field is aligned along the short (long) axis of the 
dot. We show how this can be understood from the 
form of the SO Hamiltonians and the approximate par- 
ity of the one-electron wave functions. The dependence 
of the AC width on the direction of the magnetic field 
turns to be a I sin0| shaped function and when both cou- 
plings are present this function is shifted by an amount 
which depends on the relative strength of both interac- 
tions. This shift is affected by the orientation of the 
dot within the [001] plane due to the SO bulk-induced- 
anisotropy (Dresselhaus term) . For completeness we also 
study the influence of the dot shape. We show that for a 
square-based quantum dot the anisotropic dependence of 
the AC width is only observed when both couplings are 
present — Moreover we show that for increased height of 
the dot the orbital effect of the magnetic field modifies 
the energy spectrum but the shape of the dependence of 
the anticrossing width on the direction of the in-plane 
magnetic field remains unaltered. 
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The present work is organized as follows: we start with 
an outline of our theoretical approach in section II. In sec- 
tion III we present our numerical results starting from the 
single-electron case which provides us with physical in- 
sight in the reasons for the SO coupling anisotropy. We 
continue by studying different orientations and geome- 
tries of the dot and we end the section with the two- 
electron case that allows for a direct comparison with 
the recent experimental data of Ref. [2ll We end with 
a concluding discussion in section IV and a summary in 
section V. 



II. THEORY 
A. Model 

Our aim is to calculate the energy spectra of the one- 
and two-electrons confined in a three-dimensional quan- 
tum dot in the presence of SO coupling and a magnetic 
field oriented within the quantum dot plane. The effect 
of the spin-orbit coupling on the energy are very small re- 
quiring a very high numerical precision when evaluating 
the energy spectrum. We assume that the quantum dot 
is cuboid in shape and that the confinement potential is 
separable, namely ^(r) = T4(x) -I- Vy{y) + Vz{z). More- 
over, we assume that the one-dimensional confinement 
potentials 14, Vy and Vz can be described by an infinite 
quantum well model. This is a reasonable approximation 
for not to small quantum dots. Under these assumptions 
one can construct a sufficiently precise solver for the two- 
electron problem. We consider quantum dot with varied 
in-plane orientation with respect to the crystal host. The 
z-axis is taken along the [001] crystal direction which is 
also the vertical axis of the dot. The orientation of the 
dot is described by a rotations of the x and y directions, 
which are the axes of the dot with respect to the [100] 
and [010] crystal directions. The outline of our quantum 
dot and the coordinate system used is depicted in Fig. [T] 



B. Method 

We employ the effective mass approximation with the 
single-electron Hamiltonian of the form 



z[001] 



2m* 



V{r) 



l + ^g^isB-a + HsiA + HsiA, (1) 



where k ~ — iV + 1 is the identity matrix, V{r) 
defines the confining potential and Hbia, Hsia are the 
spin-orbit coupling Hamiltonians. The x, y and z direc- 
tions are taken along the axes of the dot. But notice that 
the SO interaction Hamiltonians are defined in the coor- 
dinate system with axes parallel to the [100], [010] and 
[001] which we denote with x' , y' and z. Both coordinate 
systems are transformed into each other by an in-plane 
rotation over an angle 6. 




[010] 



[100] 



FIG. 1: Schematic of the quantum dot system with the used 
coordinate system fixed to the quantum dot. The crystallo- 
graphic directions of the InAs host lattice are also indicated 
together with the direction of the in-plane magnetic field. 



We introduce the Rashba coupling with Hamiltonian, 

Hsia = ■ {a' x k'), (2) 

where a defines the coupling strength. For infinite quan- 
tum well confinement the term V'V within the dot equals 
the external electric field. We neglect the infiuence of the 
in-plane component of the electric field24 and obtain the 
Rashba Hamiltonian in the form. 



H. 



sia 



dV_ 

a7 



ip'x' ky' (7y' kx' ) ■ 



(3) 



Thus the electric field is in the z direction which is in- 
corporated by taking a non-zero slope of the bottom of 
Vz{z). 

Inversion asymmetry of the crystal lattice results in a 
Dresselhaus SO coupling that is described by the Hamil- 
tonian 



Hbia = 7 [(^x'K'ikl - fc^O + cry,ky>{kl, - k^) 
+azkz{kl, - kl,)] , 



(4) 



where 7 is the coupling constant. 

The coordinate system used for the SO coupling can 
be transformed into the coordinate system used for the 
quantum dot through the transformation 



x' = a;cos(6') — ysin(6') 
y' = x sin{9) + y cos{9) 



(5) 



which applies both to the Pauli matrices a and the coor- 
dinates of the momentum operator. 

We include an in-plane magnetic field of orientation 
B = i? (cos 0, sin (j>, 0) which is described by the gauge 
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A = B {z ■ s'm(j),0,y ■ cos (j)). The magnetic field vector 
B for (/) = is oriented along the x direction (see Fig. [T]). 

The one-electron Hamiltonian ([T]) can be rewritten in 
the form h = + hy + + hns, where 
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2m* dx^ 



(6) 
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2m* dy'^ 



^2 a2 



2m* i9z2 



Vy{y) 



2m'' 



2m" 



■y^ cos^ (j), (7) 



-z^sin^c/), (8) 



are spin independent parts separable in the x, y and z 
directions, and 



ifieB 



;sm0— + ycos0— 



(9) 



--^dt^bB [o-o; COS (j) + Oy sin 0] + HsiA + Bbia, 



is the nonseparable part that contains the spin dependent 
terms, 

The eigenenergies and the eigenvectors 
ipx{x),ipy{y),ipziz) of the Hamiltonians hx,hy,hz 
are calculated separately on one-dimensional meshes of 
A^iD = 1000 points. In a next step we diagonalize hns in 
a basis of products of the eigenstates ip^ix)-, ipy{y), ipz{z) 
resulting in three-dimensional spin-orbitals '0(r, cr). We 
typically take = Ny = 20, Nz = 10 one-dimensional 
eigenstates (we assumed Rz <SC Rx,Ry), which including 
the degeneracy of the spin gives a basis consisting of 
8000 elements which results in an accuracy better than 
5ficV. 

We solve the two-electron problem as described by the 
Hamiltonian 



H = hi 



(10) 



using the configuration-interaction approach. In our nu- 
merical calculation we take the dielectric constant for 
InAs e ~ 14.6. Hamiltonian (jlOp is diagonalized in a ba- 
sis constructed of antisymmetrized single-electron spin- 
orbitals ^{r, <t) 



n n 



1 

* = ^E 1. [^,(l)V',(2)-^.(2)V',(l)], (11) 

where 1, 2 are the spatial (r) and spin (cr) coordinates of 
the corresponding electron. The electron-electron inter- 
action matrix element requires the calculation of integrals 
of the form. 



- (V'i(i-i)V'j (1-2)1 - 



-|^fc(i-i)'0i(r2)) 



47r£o e|ri-r2| 

; ^3 ,*/ N , / ^ / j3 e i/'*(r2)Vv(r2) 
e / d''ri?/'*(ri)^fe(ri) / (rv2- 



e / d3ri^*(ri)V^fe(ri)F,,(ri). 



- (12) 



A direct calculation of these 6 dimensional integrals re- 
quires an enormous numerical cost. Therefore, we use 
the method^^ in which the innermost integral is at- 
tributed to an electric potential Vji (ji ) originating from 
an electric charge distribution ^*(r2)^i(r2). We calcu- 
late the electric potential by solving the Poisson equa- 
tion V^Vj7(ri) = — e/(eeo)'0|(ri)^i(ri) with the bound- 
ary condition 



,3 V^*(ri)V^;(ri) 



£\rb - ri| 



(13) 



where lays within the boundary of the the computa- 
tional box. The Poisson equation is solved on a grid that 
covers the dot area. The calculation accuracy is care- 
fully monitored^^ and the configuration-interaction cal- 
culation convergence better than 10/j, eV is reached for 
n = 20. 



C. Parameters 

The bulk of our results presented in the following sec- 
tions are obtained for p arameters described below. In the 
experiment of Ref. |2lJ an anisotropic InAs sclf-organized- 
quantum-dot (SOQD) grown on a [001] GaAs substrate 
is studied with a confinement potential that is elongated 
due to the source and drain electrodes placed above the 
dot. The orientation of the dot with respect to the in- 
plane crystal directions is not well resolved and in the 
present work this is taken as an additional parameter 
which is studied. We take R^ = 100 nm as the long 
and Ry — 60 nm as the short size of the doti^i We take 
Rz = 10 nm as a reasonable estimate of the dot height 
(note that the SOQD has a nominal pyramidal shape^i 
with height 20 nm, but our model is limited to a potential 
with rectangular shape of vertical cross section). Rz in- 
fluences the effective strength of the Dresselhaus coupling 
constant and the orbital effects of the in-plane magnetic 
field. Results for i?^ > 10 nm are also discussed below. 

For the purpose of the present study it is important 
to notice that the electric field in the growth direction 
defines the strength of the Rashba coupling. The elec- 
tric field is influenced by the potential proflle within the 
dot^l the Schottky barrier at the dot/electrode inter- 
face, surface charges and applied potentials.— The elec- 
trostatics of the actual device is complex and its complete 
description is out of the scope of the present work. Never- 
theless, we are able to estimate the external electric held 
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present in the system by considering the stability dia- 
gram and the width of the systems<^ We estimated the 
maximal value of the external electric field to be of order 
—30 kV/cm for which the electrons are still present in the 
dot<^ From the gate voltage Vg = —0.4 V of two-electron 
spectroscopy we estimated Fz = —13.6 kV/cm and this 
value is used in our numerical calculation. Finally, in this 
paper we indicate that the ratio of the Rashba coupling 
strength (that is proportional to F^) to the strength of 
the Dresselhaus coupling can be extracted from the ex- 
perimentally measured orientation of the magnetic field 
for which the SO-induced AC vanishes. 

We take the SO coupling parameters as a = 1.1 nm^ 
from Ref. [13 for the Rashba coupling and 7 = 
26.9 meVnm^ from Ref. [HI for the Dresselhaus coupling 
constant. Material parameters for In As are adopted from 
Ref. [13 with values m* 0.026, g = -17.5. 

III. RESULTS 
A. Without SO coupling 
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FIG. 2: (Color online) The black solid curves present one elec- 
tron energy spectrum obtained without spin-orbit coupling for 
the in-plane magnetic field regardless of the (p value. The red 
dashed curves are the energy levels when only Dresselhaus 
coupling with 7 — 26.9 meV nm'^ is included with magnetic 
field aligned along the y direction (0 = 90°). The inset shows 
a zoom of the energy levels in the vicinity of the anticrossing. 

We consider first the dot aligned such that the x- 
axis is oriented along [100] (y-axis along [010]), namely 
9 = 0. The energy spectrum obtained in the absence 
of the SO coupling (wc take a = 7 = 0) for a single- 
electron anisotropic quantum is presented in Fig. [2] by 
the black solid curves. In the absence of the magnetic 
field the ground state is doubly degenerate with respect 
to spin and the spatial wave function is of even symmetry 
with respect to plane inversions: ip{x,y, z) = ip{—x,y,z), 
fp{x, y, z) = ^{x, -y, z), and ij{x, y, z) = ^{x, y, -z). We 
denote the state of even symmetry with respect to all in- 
versions by ['!/'+)■ The first-excited state is a spin-doublet 



with wave- functions meeting the symmetry conditions: 
tp(x,y,z) = -tlj{-x,y,z), il;{x,y,z) = ip{x,-y,z), and 
Tp{x,y,z) = %l^{x,y,—z). We will refer to this state as 
l^"-). The non-zero magnetic field lifts the spin degen- 
eracy splitting of the states of the same parity by the 
Zceman energy. The energy levels depicted by the black 
lines in Fig. [5] arc obtained regardless of the value in 
spite of the lateral anisotropy of the dot. Due to the small 
Rz value and the in-plane alignment of B, no orbital ef- 
fects of the magnetic field are observed (the influence of 
the height of the dot is studied in subsection F). 

Generally, in the presence of an in-plane magnetic 
field the Hamiltonian (1), even without SO interaction, 
does not commute with the plane inversion operators 
Px and Py [defined as Pxf{x,y,z) = f{—x,y,z) and 
Pyf(x, y, z) = f{x, —y, z)]. However, due to the insignif- 
icance of the orbital effect of the magnetic field for this 
flat quantum dot, the parity with respect to reflection 
through the x = and y = plains is approximately pre- 
served (with {Px) and {Py) above 0.97) even for non-zero 
B. For the following discussion we denote the four lowest- 
energy states for small magnetic field aligned parallel to 
the y direction as l*(+4)>' l*(-,t))' l*(-4)) 

with corresponding energies £;(+_-^) , £^(+4), £;(_^;) 
where the arrow denotes the spin state aligned parallel 
(t) or antiparallcl (J,) to the magnetic field vector. 



B. Single type of SO coupling 

Inclusion of the SO interaction lifts the spin polar- 
ization of the states and changes the crossing observed 
between the energy levels of |5'(+4)), |^'(_4)) around 
B = 4.25 T into an anti-crossing. The inset of Fig. [2] 
shows the anticrossing energy levels for = 90° (B par- 
allel to the y axis) by the red curves when only Dressel- 
haus coupling with 7 = 26.9 meVnm'^ is included. We 
denote the minimal energy difference between the anti- 
crossing levels as Eac- For applied parameters we obtain 
Eac ~ 146 /icV. Outside the anticrossing the SO interac- 
tion does not modify the energy spectrum in a noticeable 
way i.e. the black and red curves approximately coincide. 

In the presence of the SO couphng the anticrossing 
energy levels depend on the orientation of the magnetic 
field. In Figs. [3]^ a) and [3]^b) we plot the energy levels 
obtained for pure Dresselhaus and pure Rashba interac- 
tion, respectively, for three different (p values. In both 
cases clear dependence of the anticrossing width Eac is 
observed with respect to the B orientation. For pure 
Dresselhaus coupling the anticrossing is the widest when 
the magnetic field vector is perpendicular to the y direc- 
tion ((/) = 90°) [the red curve in Fig. EJa)]. When the 
field is aligned along the x direction (0 = 0) the mixing 
between levels vanishes [the black curve in Fig. [Sja)] and 
there is crossing of levels. With pure Rashba coupling 
the dependence is opposite - the anticrossing vanishes 
when B is aligned along y and Eac is largest when B is 
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FIG. 3: The energy levels near the anticrossing for pure Dres- 
selhaus (a) and pure Rashba (b) couplings for different B 
orientations. The black crosses are the results from a diago- 
nalization of the matrix (|22|) and the red crosses the results 
of diagonalization of the (|19|l matrix. The magnetic field is 
oriented along the x direction (0 — 0) for the black symbols 
and curves, and along the y direction (0 = 90°) for the red 
symbols and curves. 



aligned along x. 

The direction of the magnetic field for which the mix- 
ing between the states disappears can be infered from 
the analytic form of the SO Hamiltonians utilizing the 
approximate symmetries of the wave functions of the con- 
fined electron. Let us first inspect the case of pure Dres- 
selhaus coupling and remind that for 6 = the Hamilto- 
nian (j4|) has the same form in the x, y and z coordinate 



system. Averaging the Hamiltonian (|4]) over the z direc- 
tion one obtains, 



H 



2D 
BIA 



Oy ky] 



-7 [(Tykykl - axkxkl] 



(14) 



The second term is the so-called cubic Dresselhaus 
term which is negligible as long as the height is much 
smaller than the lateral size of the dot [i.e. until the 
value of (fc^) or (fc^) becomes comparable with (fc^)]. For 
an infinite quantum well ground-state wave function in 
the z-direction the last term in vanishes^ and 



(15) 



The simplified Dresselhaus Hamiltonian takes now the 
form 



HmA = 1^" {o-xkx - o-yky) 



..2D 



(16) 



Let us now consider the case of a magnetic field aligned 
paralel to the y direction. In our basis we include only 
the low-energy states |^(_|_ j^)), j^f (_ -!-•)) that exhibit an 
energy crossing without SO coupling. The matrix of the 
H^^js^ Hamiltonian limited to this basis is given by. 



J 



,2-D/,T/, I, ^ I. liTi. ,\ Z?. , I „,2-D/ 



l^^{'^(-,^)Wxko:-(Tyky\'^( + X)) £;(_.t) + 7 ,t) I " ^^2/ ^2/ I *(-,t) ) 

r 



(17) 



The states |^'(+ j^)), |^'(_ -j^)) are separable into an orbital on the states with definite spin one gets, 
and a spin part. Due to the action of the Pauli matrices 



{+.X)-l''°{'^{+.X)\<yyk: 
7'''(*(-,t)l^-^-l*(+,i) 



7'^(*( + ,i)|'7.fc.|*(-,t)) 

^(-,T)-7'^(*(-,t)l^ 



(18) 



r 



For the magnetic field vector aligned parallel to the y di 
rection the components of the momentum operator vec 

^,fcz = —i^z- Due to 



tor are k^ — ^~dx ~t~cBz,/cy — x-q — ^k^ — ^"Oz 



parity one obtains. 
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The non-vanishing off-diagonal matrix elements mix the 
states which results in an avoided cross- 



ing between the corresponding energy levels. By the red 
crosses in Fig. [3] we plot numerically calculated eigenval- 
ues of the matrix ([T^. Note that the crosses and lines 
are in perfect agreement proving that for our dot with 
the assumed geometry the H^j^ is in fact a good ap- 
proximation to Hbia- 
dx\'^ (-A)) \ Let us now consider the case of a magnetic field aligned 
/parallel to the x (</> = 0) direction. In this case the 
(19) low-energy states which energy levels cross without SO 
coupling are |\I'(_|_^)), |^(_^^)), where the arrow denotes 
the electron spin aligned parallel (— >) and antiparallel 
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(<— ) to the magnetic field vector B. The matrix of the 



Due to spin 



For the magnetic field aligned along the x direction the 
components of the momentum operator vector are, = 
~^m^^v = = + ^^y- integrals in Eq. 

(PT|) vanish due to the parity of the states and we finally 
obtain, 



The matrix consists only of diagonal elements that 
are equal to the energy of the basis states. Thus the 
|^(_^^)) states are not mixed by the Dressel- 
haus coupling in this configuration and there is no anti- 
crossing of energy levels. We plot the eigenvalues of the 
matrix ([22) by the black crosses in Fig. [3{a). 

A similar analysis can be made for the Rashba Hamil- 
tonian Due to the fact that the analytic form of both 
Hamiltonians Hsia and Hg^j^ are similar, i.e. only the 
kx and ky are swapped (and the coupling constants are 
different), it is clear that the dependence of AC width on 
magnetic field direction is opposite - the mixing between 
the states vanishes when the magnetic filed is aligned 
along the y direction. 

C. Anisotropy in the presence of both SO couplings 

Let us now consider the effect of both Dresselhaus and 
Rashba coupling. Figure |4] presents the avoided cross- 
ing energy Eac a-s a function of the angle (j) between 
the X axis and the magnetic field. For pure Dressel- 
haus (the blue dashed curve in Fig. |4]) and pure Rashba 
(the green dotted curve in Fig. |4]) coupling the extrema 
are shifted by 90° in agreement with our previous anal- 
ysis. The curves in Fig. |4] are accurately described by 
I sin((/) — (t>Ac)\ which is the same functional form as the 
one observed in the experimental work of Ref. [2l| in Fig. 
3(f) [where the behavior was described by | cos{(j)— 0o)|]. 



HajA Hamiltonian in this two state basis is 



(20) 



(21) 

I 

200 1 ^ 1 ^ 1 ^ 1 ^ 1 




0^^^ 90 180 270 360 
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FIG. 4: The spin-orbit-induced anticrossing width Eac for 
pure Dresselhaus (blue dashed curve), pure Rashba (green 
dotted curve) and both (red solid line) the interactions 
present. For 7 = the magnetic field is -B = 4.268 T, for 
the two other cases B = 4.277 T. 

Moreover the maximal value of Eac is of the same order 
as the magnitude observed experimentally. When both 
SO interactions are present the dependence of the an- 
ticrossing width is plotted in Fig. S] by the red curve. 
The shape of the latter is the same as for pure Dressel- 
haus/Rashba coupling with pronounced minima where 
Eac is zero. When the magnetic field is aligned along 
the X or y direction the Eac equals the value for pure 
SO coupling. Note that the maxima are larger than the 
ones observed for pure couplings and its positions are 
now shifted and are no longer aligned along the axes of 
the dot. For a = 1.1 nm^ and 7 = 26.9 meVnm'^ the 
shift of the dependence is cf>AC = 27.8°. The latter value 
can be understood as follows. Let us denote the direc- 
tion of the magnetic field for which the AC vanishes for 
pure Dresselhaus and pure Rashba couplings by the vec- 
tors dsiA and dsiA, respectively. Next, we estimate 
the strength of each interaction. Maximal induced anti- 



l^°{'f{-^^)Wxkx - cryky\^(+,^)) + l^^{'^{-.^)Wxkx - cryky\^(_ 



one gets. 
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crossing width is E'^y^ = 146 /xeV and E%1^ = 77 ^eV 
for Dresselhaus and Rashba coupling, respectively. Thus 
the Dresselhaus interaction is 1.9 times larger than the 
Rashba coupling what makes the vector Abia 1-9 times 
longer than dsiA- Let us denote the magnetic field for 
which the effect of both spin-orbit couplings is zero by 
the vector Abia+sia = Abia + ds/A- It is easy to show 
that this vector forms an angle = 27.8° with the x axis. 
Thus when both couplings are present, the effect of the 
total spin-orbit coupling disappears when the external 
magnetic field is directed along this vector. In fact that 
is exactly what we observe in our calculation (see posi- 
tion of the minimum of the dependence depicted with 
the red curve in Fig. |4]). The formula | sin(0 — 0ac)| 
refiect the fact that the dependency obtained for both 
SO couplings present can be considered as an absolute 
value of a sum of the dependencies obtained for pure SO 
couplings described by — cos 4> and sin (p for pure Rashba 
and Dresselhaus couplings respectively. 



D. Dependence on the quantum dot orientation 
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FIG. 5: Avoided crossing energy width as function of the di- 
rection {(j)) for different orientations of the dot with 6 = 10° 
(black curves), 9 = 30° (blue curves), 9 = 45° (green curves), 
9 — 60° (violet curves), 9 = 80° (red curves). Results are 
shown for (a) pure Rashba, pure (b) Dresselhaus and (c) for 
both the couplings present, (d) The value of the magnetic 
field angle ^ac at which Eac = as a function of the an- 
gle 9 for pure Rashba (green diamonds), pure Dresselhaus 
(black squares) and both couplings present (red dots). The 
red curves are obtained from Eq. H23p . 

Different in-plane orientations of the anisotropic po- 
tential of the dot with respect to the crystal host arc 
now considered where the long axis of the dot forms an 
angle 6 with [100]. In Figs. [Sl[a,b,c) we present the size of 
the avoided-crossing as a function of the direction of the 
rotated magnetic field (note that the angle is defined 
as an angle between the magnetic field vector and the 
long axis of the dot) for six different orientations of the 
dot. The dotted curves in Fig. [5ja) presents the result 
obtained for pure Rashba coupling. We observe that the 
Eac dependencies are exactly the same as in Fig. U re- 
gardless of the dot alignment. The minimum of the Eac 



does not change its position and the energy levels are not 
affected by the orientation of the dot. We show in Fig. 
[5jd) the <j>AC angle for which the Eac = as a function 
of the angle 9 by the green diamonds. 

For pure Dresselhaus coupling the Eac dependencies 
[depicted by dashed curves in Fig. EKb)] arc shifted as the 
dot is rotated. For the case studied in previous subsec- 
tions (where 9 = 0) the AC vanished when the magnetic 
field was aligned along the long axis of the dot (0ac = 0). 
When the dot is oriented by 9 = 45° (long axis oriented 
along the [110] direction), the anticrossing vanishes when 
the magnetic field is aligned along the short axis of the 
dot [see green dashed curve in Fig. ^h)] ~ (pAC = 90° - 
the same as for pure Rashba case. We plot in Fig. [S^d) 
the angle 4>ac for pure Dresselhaus coupling by the black 
squares for different orientations of the dot. We find that 
the angle exhibits a (j>AC = —29 dependence [black solid 
lines in Fig. O^d)]. Moreover we observe that for both 
cases when only a single type of SO coupling is present 
the maximal value of the AC width remains unchanged. 




[100] 



FIG. 6: Schematic of the method of calculation of the angle 
0AC for which the the AC vanishes when both SO couplings 
are present and the quantum dot (dashed rectangle) is ori- 
ented with its long axis forming the angle 9 with the [100] 
direction. The vectors depict the directions of B for which 
Eac = for pure Rashba coupling (green arrow), pure Dres- 
selhaus coupling (blue arrow) and both couplings present (red 
arrow). The coordinate system connected with the dot axes 
and the crystallographic directions is also shown. 

In Fig. [5{c) we show the results when both SO cou- 
plings are present by the solid curves. The maximal val- 
ues of Eac and the angle 4>ac for which the minima 
are observed change when the dot orientation is varied. 
Both facts can be understood similarly as discussed in 
subsection C. We can justify the (j)Ac values considering 
the orientation of the (Ibia+sia = Abia + Asia vector. 
But now the orientation of the Abia vector assigned with 
Dresselhaus coupling is changed as the dot is rotated, ie. 
the vector Abia forms an angle —29 with the long axis 
of the dot. The rotation of the dot does not change the 
maximal value of Eac when only a single type of SO cou- 
pling is present and the previously derived value for the 
relative strength of both couplings remains unchanged 
(and thus also the ratio of the length of the Abia and 
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d-siA vectors). We take 1 as the length of dsiA and 1.9 
as the length of dsiA- In Fig. [S]we schematically present 
the considered vectors and the angles they form with the 
axes of the dot. The angle between the dBiA+siA vector 
(red arrow in Fig. |6]) and the x direction can be easily 
calculated 



(pAc = arctan 



1 + 1.9sin(-26') 
1.9cos(26') 



(23) 



With the red dots in Fig. [5Kd) we plot the angle (pAC 
obtained from our numerical calculation in the presence 
of both couplings for different orientations of the dot 
which agree very well with the values (red curves) ob- 
tained from Eq. (|23p . Along with the changes of the 
orientation the length of the dsiA+siA vector is changed 
which results in different values of the maximal AC width 
observed in Fig. EJc). 

A systematic study of the value of the (pAc angle de- 
pendence on the SO coupling strengths and the dot align- 
ment is given in subsection G where the two-electron case 
is studied. 




FIG. 7: Energy levels of one-electron quantum dot with 
square base with Rx — Ry — 100 nm. The black curves cor- 
respond to = 0, the blue curves to (f> = 45°, yellow dotted 
to (j) — 22.5° and the red dashed curves to (/> = 90°. In (c,d) 
we additionally plot the energy levels obtained for p — 135° 
with green dotted curves, (a) Pure Dresselhaus coupling, (b) 
pure Rashba interaction, (c) both SO interactions are present 
and (d) both SO interactions are present with a increased by 
a factor of two. 



F. Larger dot height 



E. Quantum dot with square base 

The above discussion was for a lateral anisotropic 
quantum dot. Now we study the case of a dot with sym- 
metrical base (we assume Rx = Ry = WO nm) and 9 = 
and investigate if this has an influence on the anisotropy 
induced by the SO coupling. In the absence of the SO 
interaction and a magnetic field the first-excited state is 
spin-doubly degenerate due to parity. The magnetic field 
lifts the spin degeneracy but the degeneracy due to par- 
ity is not removed. The inclusion of a single type of SO 
interaction induces a repulsion between the energy levels 
of the ground-state and one of the states from the parity 
doublet [see the red dashed curves in Fig. [TKa) for the 
case of pure Dresselhaus coupling and Fig. [7Kb) for pure 
Rashba coupling] . The same configuration of energy lev- 
els is obtained regardless of the angle 0. In both Figs. 
[TJa) and Ulh) the black (</> = 0), blue (0 = 45°), yel- 
low dotted {(j) = 22.5°) and red dashed curves {(f> ^ 90°) 
coincide. The dependence of the energy levels on the 
magnetic field orientation starts to appear already when 
the dot is elongated by a factor of 1%. 

However, when both Rashba and Dresselhaus interac- 
tions arc present the AC width varies with the rotation 
of the magnetic field - Figs. [7| (c,d). We observe that 
the anisotropy is most pronounced when a is increased 
by a factor of two - the case when both couplings have 
comparable strengths i^i^ In such a case when the mag- 
netic field is directed along the diagonal, i.e. (j) = 45° [the 
blue curves in Fig. [T^d)] the anticrossing between the en- 
ergy levels of the ground-state and both states from the 
parity-doubled becomes very small. 
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FIG. 8: One-electron energy levels for rectangular-based 
quantum dot with Rz = 20 nm (a) and Rz = 40 nm (b) 
for (f) = 0, (j) = 45° and <j) ~ 90° plotted with the black solid, 
blue solid and red dashed curves, respectively. The red dots in 
the insets to both plots presents the anticrossing width Eac 
obtained from the energy-spectrum for a given (j> value and 
the black curves are the fitted | sin{(j> ~ 4'ac)\ dependencies. 

Let us now return to the case of the quantum dot with 
rectangular base with R^ = 100 nm and Ry = 60 nm. 
For the previous dot with R^ = 10 nm no orbital effects 
from the magnetic field on the energy spectrum was ob- 
served [see the black curves in Fig. [2] . However this is no 
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longer true for larger Rz values. This can be seen from 
Figs. [Hl^a) and |H{b) where we plot the energy levels of a 
quantum dot with height = 20 nm and R^ — 40 nni, 
respectively, in the presence of SO coupling (with both 
SO interactions present). The energy levels depend on 
the magnetic field orientation even outside the anticross- 
ing region. This is due to the elongation of the confine- 
ment potential in the x direction. The SO-induced an- 
ticrossing is shifted to lower magnetic fields as the value 
of the angle cj) becomes closer to (j> = 90° [this is anal- 
ogous to the experimental observation - compare with 
Fig. S7(a) from Ref. [2^. We calculated the anticrossing 
widths for different values of 4> and plot them as red dots 
in the insets of Fig. [8l^a,b). Then we fitted the points 
with the function A\ sin(0 — 0ac)| where A = 86 /xeV, 
= 65° for Rz = 20 nm and A = 77 iieV, (t)AC = 82° 
for Rz ~ 40 nm. Notice the agreement between the fitted 
curve and the data points. From this fact we conclude 
that in spite of the presence of orbital effects, previously 
found dependence of the anticrossing width on the angle 
(/) still holds, but with modified A and 4>AC values. The 
latter fact can be attributed to the reduction of the Dres- 
selhaus coupling strength. This can be accounted for by 
considering the Dresselhaus coupling Hamiltonian p6)) 
in which the coupling strength decreases as {1/Rz)'^. In 
the calculation performed for pure Dresselhaus interac- 
tion we obtain the maximal Eac values 146 /ieV, 37 /leV 
and 11 AieV for Rz = 10 nm, Rz = 20 nm and Rz = 40 
nm, respectively. For increased dot height the obtained 
Eac values decrease approximately as (1/i?^)^ with the 
largest discrepancy for large Rz value (i.e. when the ap- 
proximation of the coupling strength by Eq. be- 
comes inaccurate) . The decrease of the Dresselhaus cou- 
pling strength for increased height of the dot results in 
a shift of the Eac dependency on </> towards the one 
obtained for a fiat quantum-dot with only Rashba in- 
teraction present (compare the black curve in the inset 
of Fig. El^b) with the green dotted curve in Fig. |4]) - 
(/)Ac becomes close to 90°. Also the maximal Eac value 
becomes closer to the one obtained for pure Rashba cou- 
pling - A tends to 77 /ieV with increasing Rz- The shift 
in the ^ac value [see insets of Figs. [5] (a,b)] can be un- 
derstood from the relative strengths of the Rashba and 
Dresselhaus coupling as discussed in section III. C. 



G. Two electron results 

In a recent experiment [2l| the ground-state and ex- 
cited states were measured provided that the latter en- 
tered into a finite but narrow transport window deter- 
mined by the voltages applied to the source and drain 
electrodes. The avoided crossings that appear for a 
single-electron in the excited part of the spectrum, which 
we described above, were outside the transport window. 

In the two-electron regime and in the absence of both 
the magnetic field and the SO interaction the ground- 
state is a spin-singlet and the first excited state is a spin- 



triplet. Under the presence of an external magnetic field 
the ground-state singlet energy crosses the triplet energy. 
When we turn on the SO coupling it induces an avoided 
crossing between the states of opposite spin which was 
well resolved in the experiment |2l| . 



310 I ' ^ I I 




12 3 
B[T] 

FIG. 9: Two-electron energy spectrum in the presence of both 
Rashba and Dresselhaus SO coupling, for angles </> = (black 
curves), (\> = 27.8° (green curves), = 45° (blue curves) and 
= 90° (red curves). The inset shows the energy levels in 
the vicinity of the anticrossing. The results are obtained for 
61 = 0. 

Similarly to the one-electron case the SO coupling is 
responsible for changes in the size of the anticrossing en- 
ergy when the orientation of the magnetic field is varied. 
Figure [3] presents the low-energy spectrum of the two- 
electron quantum dot in the presence of both Rashba 
and Dresselhaus coupling for a dot aligned with its long 
axis along the [100] direction (d = 0). In the inset we 
plot the energy levels in the vicinity of the avoided cross- 
ing. The anticrossing vanishes for exactly the same angle 
(I) AC = 27.8° as for the one-electron case discussed above 
(see the green curves in the inset of Fig. 

In Fig. [Tola') we plot the angular dependence of the an- 
ticrossing width Eac for pure Dresselhaus, pure Rashba 
and when both couplings are present by the blue dashed, 
green dotted and solid red curves, respectively. Notice 
that all three dependencies have the same shape as for 
the case of the one-electron considered in subsection B 
(compare with Fig. 2]), only the maximal Eac values are 
about 1.5 times smaller. 

As was presented in subsection C for the anisotropic 
quantum dot the angle (pAC depends on the relative 
strength of both SO interactions and the in-plane orienta- 
tion of the dot (explained in subsection D). On the other 
hand the 4>ac value can be measured experimentallji^i 
and the orientation of the quantum dot with respect to 
the crystal directions can be obtained by inspecting the 
facets of the dot. This opens the possibility to employ 
such a measurement to evaluate the relative strength of 
the Rashba and Dresselhaus couplings for a dot with 
given orientation with respect to the crystal host. Let 
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FIG. 10: (a) Width of the singlet-triplet avoided crossing as 
a function of the angle (j) for pure Dresselhaus coupling (blue 
dashed curve), pure Rashba coupling (green dotted curve) 
and for both couplings present (red solid curve). The inset 
shows the experiment results (symbols) of Ref. [SJ together 
with the result of the present calculation (purple curve) with 
7 = 29.58 meVnm'' and a — 4.731 nm^. The magnetic field 
is B = 2.211 T for pure Rashba coupling and B = 2.209 T 
for all the other cases. The dot is aligned with its long axis 
along [100] i.e. {6 = 0). (b) The 0ac value for different 
strength ratios of the Rashba and Dresselhaus coupling for 
four different orientations of the dot. The symbols present the 
results of our numerical calculation and the curves represent 
the analytical result given by Eq. (|24|l . The circles show the 



results obtained for Ry 



60 nm, Rz 



10 nm, the crosses 



10 nm, the triangles for Ry = 20 nm, 
Rz = 10 nm, and the diamonds for Ry — 60 nm, Rz = 20 nm. 
In all cases Rx is 100 nm. 



us define the strength ratio of the SO interactions as the 
ratio of tlie effective coupling constants a* and . The 
Rashba coupling strength denoted with a* is calculated 
as a* = —a [^] = —a\e\Fz and the Dresselhaus cou- 
pling 7^^ is obtained from Eq. (|15p . 

We previously derived the angle (j)Ac for given relative 
strength of the SO couplings for a given orientation of 
the dot [see Eq. (P^ j. Let us substitute the 1/1.9 value 
by a*/7^'^ in Eq. ((231) from which we obtain 



a 

^2D 



a_\e\FJ^ 

7 TT-^ 



cos(26l) [tan( 



tan(-20)] . 

(24) 



This function is shown in Fig. [TUf b) by the solid lines 



for different orientations of the quantum dot. With the 
black symbols we mark the angle ipAC obtained from our 
numerical calculations for dots with different geometries 
(see figure caption) with = for different SO coupling 
strengths. For such case with 6* = (the dot oriented 
with its long axis along [100]) and pure Dresselhaus cou- 
pling (a*/7^^ = 0) we obtain (j)AC = 0. When the 
Rashba coupling strength is increased the points move to- 
ward the angle (j)AC = 90° obtained for pure Rashba SO 
coupling. The green, red and blue symbols in Fig. rTOf b) 
are the (jiAC values obtained from our two-electron nu- 
merical calculation for different orientation of the quan- 
tum dot. 

In the above discussion we assume that a* — ~a\e\Fz 
and 7^^ — -fTr'^/Rl describe the strength of the spin- 
orbit interactions. For the Rashba coupling given by 
the Hamiltonian ([3]) [i.e. when an electric field is only 
present in the growth direction] the above a* expression 
is valid regardless of the dot geometry. However, due to 
the fact that 7^^ originates from the Hamiltonian 
it describes the strength of the Dresselhaus coupling cor- 
rectly only when the cubic term 7 ^cTykyk^ — axkxky^ is 
negligible, which is the case when , Ry ;» R^ and when 
the term with (kz) is close to zero i.e. for a dot with 
limited height 1^ All the symbols in Fig. [TUT b) approxi- 
mately coincide with the dependency given by Eq. (|24p . 
Discrepancy is seen in the limit of a narrow dot with 
Ry = 20 nm (the triangles) and for increased height of 
the dot for R^ = 20 nm (diamond symbols). We conclude 
that for anisotropic quantum dots with limited height the 
Q!*/7^-° ratio is a good measure of the relative strength 
of the Rashba and the Dresselhaus spin-orbit couplings 
which can be estimated from analytic expression (|24p . 

The experiment of Ref. [2l| found (j)Ac = 59° and we 
can use Eq. to calculate the relative strength of 

the SO interactions. However, as the orientation of the 
anisotropic potential of the dot with respect to the crys- 
tal directions was not resolved in the experiment we need 
to assume a value for 9. We take 6 = and by matching 
the absolute value of the SO coupling constants (through 
the maximal value of Eac) we obtained a* /-f"^^ ~ 1.66 
by fitting the experimentally measured values for the AC 
width with our simulation results. In the inset to Fig. 
[TUfa) we plot our results (purple curve) for the SO cou- 
pling constants 7 = 29.58 meVnm'^, a = 4.731 nm^ to- 
gether with the data points from Ref. [2l|. However as the 
relation between the crystal directions and the long axis 
of the dot is not known the fit only proves the validity of 
the discussed process behind the anisotropy and not the 
exact value of the ratio a*/7^^. Moreover as the electro- 
statics of the actual device is complex the presented result 
is not the exact simulation of the experiment. Therefore, 
we present in Table I the strength ratios for different ori- 
entations of the dot. Note that Eq. ([M)) does not allow to 
calculate the relative strength of the couplings for a dot 
aligned with long axis exactly along [110] or [110]. In a 
such configuration for pure Dresselhaus as well as for pure 
Rashba coupling the AC vanishes for (j)Ac = 90° [com- 
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1.66 


40° 


1.27 


45° 




75° 


-0.94 


90° 


-1.66 



TABLE I: Calculated strength ratios of the SO couplings for 
<j>AC ~ 59° and different orientations of the dot. 

pare dotted curves in Fig. ^a.) with green-dotted curve 
in Fig. [5Jb)] and by that for both couplings present si- 
multaneously the minimum of the Eac dependence on (jj 
is not shifted irrespective of the coupling strength ratio. 



IV. DISCUSSION 

In the present paper we discussed the avoided-crossings 
of energy levels as induced by the presence of differ- 
ent SO coupHngs. Only for the case of a square-based 
quantum-dot [see Figs. [7Ja,b)] the dependence of AC 
width as function of the magnetic field direction was 
observed solely for both couplings present with compa- 
rable strength. This result is related to those of Ref. 
[23I where the spin-splitting of sm^/e-electron energy lev- 
els in strictly iu;o-dimensional circular quantum dots in 
the presence of a small in-plane magnetic field (before 
the crossings/ avoided-crossings appear) was calculated. 
When Drcsselhaus and Rashba coupling strengths are 
equal a well known high symmetric case is found which is 
beneficial for many spintronics applications 1^*^^ For that 
special case the energy spectrum is not aflFected by SO 
interaction effects and the spin in the [110] direction is 
strictly defined. The Zecman interaction lifts this sym- 
metry and results in a spectrum that depends on the 
orientation of the magnetic field as discussed in Ref. [l^. 
Since for equal coupling strengths the spins in the [110] 
direction are well defined, the Zeeman interaction for B 
oriented along [110] does not produce any AC between 
energy levels of spin-orthogonal states [see the blue curve 
in Fig. nid)]- 

On the other hand, in the presence of a vertically ori- 
ented magnetic field, the size of the Zeeman interaction 
induced lifting of the symmetry depends on the in-plane 
orientation^ and also the widtli^ of the dot what re- 
sults in changes in both the AC width and the effective 
g-factor which are solely observed when both SO inter- 
actions are present with comparable strength. However, 
changing the dot orientation is hardly achievable experi- 
mentally and therefore in the present work we considered 
an anisotropy that can be probed by changing the orien- 
tation of the magnetic field. 

In the present work we investigated the anisotropic de- 
pendence of the avoided-crossing width that occurs even 
for a single type of SO coupling [sec Figs. [HJa) and|3Jb)]. 
This effect is strictly connected both with the elongation 



of the confinement potential and the in-plane alignment 
of the magnetic field [see the discussion in subsection B] . 
The exact shape of the confinement potential is not im- 
portant for the studied phenomena which is a generic 
propriety of a spin-orbit-couplcd quantum dot. In our 
analysis we indicated the trends that determined the de- 
pendence of Eac on cj), in particular the dependence on 
the dot geometry [for the dot with increased height and 
for different lateral sizes of the dot - see the black sym- 
bols in Fig. [TU]that in spite of the different geometries of 
the dot still undergo the same analytical dependence Eq. 
([M)) ] or the orientation of the quantum dot with respect 
to the crystallographic directions (which infiuences the 
position of the minima of Eac purely due to Dresselhaus 
coupling - see discussion in section III. D). 

The present study shows that for an elongated quan- 
tum dot with pure Rashba coupling the anticrossing van- 
ishes always when the magnetic field is aligned along the 
short axis of the dot [see the minima of the dotted curves 
in Figs, mma) andfTUra)]. Only the presence of Dressel- 
haus coupling can result in a (pAc value that is different 
from 90°. The magnetic field direction {(pAC = 59°) for 
which the anticrossing vanished in the experiment of Ref. 
I21I suggests both SO couphngs are present, contrary to 
the argumentation provided in Ref. [2l|. The authors 
suggested that the Drcsselhaus coupling would not in- 
duce mixing between the two lowest-energy states due 
to their well defined and different values of the total an- 
gular momentum J_ = L — S* in a high magnetic field. 
However we found, that due to the in-plane alignment of 
the magnetic field^ the Dresselhaus coupling in fact in- 
duces avoided-crossings in the energy spectrum of a flat 
quantum dot [see Fig. [5] and leads also to a shift in 
the dependence of the AC width on the magnetic field 
direction [see Fig. O Fig. [5l^a,b) and insets to Fig. [5]. 



V. SUMMARY AND CONCLUSIONS 

We presented a study of the energy spectrum of one 
and two-electron spin-orbit-coupled three-dimensional 
quantum dots in the presence of an external in-plane 
magnetic field. We found that the size of the avoided- 
crossings in one- and two-electron energy-spectrum oscil- 
lates as a function of the orientation of the magnetic field. 
The oscillatory behavior could accurately be described 
by I sm{<j) — (l)Ac)\ which agrees with recent excited-state 
spectroscopy measurements performed on InAs gated 
self-organized-quantum doti^ 

For a quantum dot which is elongated in the [100] 
direction and when only a single type of SO coupling 
is present the avoided crossing vanishes for 4)ac ~ 
(0AC = 90°), i.e. when the magnetic field is aligned par- 
allel to the long (short) axis of the dot for Dresselhaus 
(Rashba) coupling. We explain this behavior as a conse- 
quence of parity- and spin-dependent mixing of the states 
caused by the SO interaction. When both couplings are 
present the (pAC value varies between and 90° and the 
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ratio of the relative strength of the interactions follows 
a tan{(pAc) dependence. The change of the in-plane dot 
orientation results in a change of (j)AC which is observed 
only when Dresselhaus coupling is present. We show that 
the experimentally measured (/)ac valued along with the 
knowledge of the orientation of the dot can be used to 
determine the ratio of the strengths of the individual SO 
interactions in case of anisotropic quantum dots. 
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